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ABSTRACT 

Piezohydraulic  actuation  is  the  use  of  hydraulic  fluid  to  rectify  the  high  frequency,  small  stroke  oscillation 
movement  of  the  piezoceramic  stack  into  a  unidirectional  movement.  This  paper  presents  a  quasi-steady  incom¬ 
pressible  flow  model  and  a  compressible  flow  model  for  performance  simulation.  These  two  models  are  derived 
based  on  the  energy  equations  for  hydraulic  flow  in  circular  pipe.  Major  losses  and  minor  losses  are  identified 
and  incorporated  in  the  model.  Assumptions  and  approximations  were  made  to  minimize  computation  effort 
while  achieving  good  accuracy.  A  piezohydraulic  actuation  system  with  active  valves  is  then  built  and  tested. 
Comparison  with  test  results  shows  that  the  simulations  accurately  predict  system  performance  under  loads  at 
which  value  leakage  is  not  present.  The  results  reveal  that  friction  losses  due  to  viscosity  are  a  major  limitation 
of  performance  for  the  current  test  setup  when  operating  at  higher  frequencies. 

1.  INTRODUCTION 

Piezoelectric  ceramics  can  produce  large  stress  but  are  limited  in  their  strain  output.  Conventional  materials 
can  produce  blocked  stresses  on  the  order  of  25-50  MPa  at  field  levels  of  1-2  MV/nr,  but  the  strain  output  is 
limited  to  0.1%  to  0.15%.  Newer  single-crystal  ceramics  produce  larger  strain  levels  (  1%)  but  lower  stresses 
(  10-20  MPa)  at  equivalent  field  levels. 

For  the  last  two  to  three  years  several  groups  have  studied  piezohydraulic  actuation  as  a  means  of  overcoming 
the  strain  limitations  in  piezoelectric  ceramics  without  sacrificing  force  capability.  The  first  use  of  piezoelectric 
materials  for  producing  1-10  W  mechanical  power  output  was  reported  by  Mauck  and  Lynch  (1;  2).  They  discuss 
the  development  of  a  system  consisting  of  a  reciprocating  piezoelectric  pump  and  passive  check  valves.  They 
report  a  no-load  speed  of  72.5  nnn/sec  and  a  blocked  force  of  271  N  for  a  particular  value  of  cylinder  bore 
diameter  (3).  Assuming  a  linear  relationship  between  these  parameters,  the  maximum  mechanical  power  output 
is  estimated  to  be  5  Watts.  More  recently,  Regelbrugge  and  Anderson  numerically  analyzed  the  performance  of 
a  piezohydraulic  system  and  highlighted  the  need  for  systems- level  analysis  (4) .  Further  experimental  work  was 
performed  by  Sirohi  and  Chopra  (5).  A  maximum  no-load  velocity  of  1.2in/sec.  and  block  force  of  351bs  was 
reported. 

Our  work  in  piezohydraulic  actuation  has  focused  on  the  development  of  systems  that  incorporate  active 
valves  for  flow  control.  Our  previous  modeling  efforts  have  concentrated  on  the  development  of  lumped  parameter 
models  of  piezohydraulic  systems  (6)  and  experimental  demonstration  of  a  closed  hydraulic  circuit  (7;  8).  Our 
experimental  results  demonstrate  the  importance  of  valve  timing  and  fluid  compressibility  in  the  optimization 
of  the  output  power  of  the  system.  These  results  qualitatively  match  our  numerical  analysis  of  the  effects  of 
compressibility  on  the  efficiency  of  piezohydraulic  and  piezopneumatic  systems  (9). 

In  this  work  we  present  quasi-static  steady  flow  models  based  on  different  approximations  of  the  hydraulic 
fluid.  In  the  first  model,  we  assume  a  fully-developed  incompressible  viscous  flow  and  in  the  second  model  we 
incorporate  the  compressibility  of  hydraulic  fluid  into  our  model.  These  two  models  are  derived  based  on  the 
energy  equations  for  hydraulic  flow  in  circular  pipe.  Major  losses  and  minor  losses  are  identified  and  incorporated 
in  the  model.  Assumptions  and  approximations  were  made  to  minimize  computation  effort  while  achieving  good 
accuracy.  A  prototype  piezohydraulic  system  with  active  valves  was  then  built  and  tested. 
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Figure  1:  Bench  Top  Experimental  Test  Setup 
Table  1:  Status  of  stack  and  valves  at  four  stages. 
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2.  SYSTEM  MODELS 

A  typical  piezohydraulic  pump  consists  of  five  major  components:  a  stack  actuator,  a  pumping  chamber,  an 
accumulator,  a  hydraulic  cylinder  and  valves  and  fittings.  Figure  1  shows  the  prototype  of  the  piezohydraulic 
pump.  The  active  valves  are  used  to  control  the  flow  of  the  hydraulic  fluid.  An  accumulator  is  used  to  provide 
the  system  with  a  bias  pressure. 

The  model  of  the  system  will  be  introduced  based  on  the  four  stage  operation  fashion  (9).  The  stack  actuator 
moves  in  a  reciprocal  fashion  under  periodic  electrical  excitation.  The  operation  of  the  system  can  be  divided 
into  four  stages:  compression,  exhaustion,  expansion  and  intake  based  on  the  phases  of  valve  operation  related 
to  the  stack.  Table  1  shows  the  status  of  valves  and  stack  in  each  stage.  For  this  four  stage  model,  it  is  noticed 
that  only  in  exhaustion  stage  does  the  system  perform  work  on  the  load. 

Laminar  flows  and  turbulent  flows  are  very  different  flow  patterns  which  are  determined  by  the  Reynolds’s 
number  of  the  flow.  For  flow  in  pipes,  it  is  laminar  when  Reynolds’s  number  is  less  than  2000  and  flow  is  usually 
turbulent  when  it  is  larger  than  3000.  The  critical  operation  frequency  for  laminar  flow  was  determined  in  (10) 
to  be  fcr  =  960  Hz.  For  operation  frequency  below  960  Hz,  the  flow  can  be  considered  as  laminar. 

2.1.  Incompressible  viscous  flow  (IVF)  model 

When  developing  this  flow  model,  we  assume  an  incompressible  flow  in  the  open  chamber  operations:  the 
exhaustion  stage  and  intake  stage,  while  the  compressibility  is  considered  in  the  closed  chamber  operations:  the 
compression  stage  and  expansion  stage. 

The  electromechanical  coupling  model  of  the  piezoceramic  stack  will  be  incorporated  into  the  flow  equations, 
which  enable  us  to  simulate  response  under  a  certain  characteristics  of  a  power  output.  Nasser  and  Leo  (9) 
presented  a  one  dimensional  model  for  piezoceramic  stack  with  many  layers  of  piezoelectric  elements: 

x  =  x0V-  -^F,  (1) 

a 


Q  =  CV-  x0  F. 


(2) 


(a)  Control  volumes  in  the  exhaustion  stage 


(b)  Control  volumes  in  the  intake  stage 


Figure  2:  Control  volumes  for  IVF  model 


where  x  is  the  displacement  of  stack,  V  denotes  the  voltage  applied  to  the  stack,  Q  represents  the  charge  and  F 
is  the  force  generated  by  the  stack,  which  is  equal  to  the  pressure  times  with  cross  section  area  of  chamber  Ach- 
The  variable  x-o  denotes  the  displacement-to- voltage  constant  of  the  stack,  Ka  is  the  stiffness  and  C  is  stress-free 
capacitance  of  the  piezoceramic  stack.  This  coupling  model  will  be  incorporated  into  the  flow  model  which  will 
be  developed  in  the  following  sections. 

At  beginning  of  a  cycle,  we  assume  that  the  system  is  pre-pressurized  to  Pi ,  the  initial  volume  in  the  pumping 
chamber  is  Vo  and  the  displacement  of  the  piston  is  xi.  Then  the  system  begins  to  operate  in  the  four  stage 
fashion.  The  models  developed  for  each  stage  can  be  integrated  to  simulate  the  full  cycle  operation  of  the  system. 

Compression  Stage 

In  this  analysis,  we  assume  that  the  bulk  modulus  of  the  fluid  j3  is  a  constant,  which  is  defined  by 


/9  =  - 


dP 

dV/V 


(3) 


Integrating  the  above  equation  and  substituting  in  the  initial  conditions  yields 

Vp  _ 

P  =  P1^/3  1n^f - ,  (4) 

j&~Xl 


Exhaustion  Stage 

We  assume  the  initial  states  of  the  stage  is  x  =  X2  and  P  =  P2.  In  this  stage,  the  fluid  flows  from  the 
pumping  chamber  into  the  high  side  of  the  cylinder.  And  the  piston  in  the  cylinder  pushes  the  fluid  to  flow  into 
the  accumulator. 

The  energy  equation  of  fully-developed  steady  laminar  flow  in  a  pipe  is  derived  by  White  (11) 
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where  subscript  i  and  2  denotes  the  physical  variables  at  different  sections  of  the  pipe.  For  a  laminar  flow  the 
coefficients  of  velocity  profile  is  a  constant:  oti  =  a2  =  2.  To  apply  the  above  energy  equation,  control  volumes 
are  chosen  as  the  dash  lines  in  Figure  2(a). 

To  simplify  the  computation,  following  geometry  constants  are  defined.  A  is  defined  as  the  area  of  chamber 
to  the  piping  and  7  is  defined  as  area  ratio  of  the  hydraulic  cylinder.  For  a  single-ended  cylinder  7  is  larger  than 
1  and  for  a  double-ended  one  it  is  equal  to  1.  The  variable  e  is  defined  as  the  area  ratio  of  the  hydraulic  cylinder 
to  the  pipe.  Substitute  the  above  constants  into  Equation  (5)  and  simplify  to  yield 


e2KL  1 
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2 


(6) 


In  the  same  manner,  we  can  derive  the  flow  equations  for  the  control  volume  CV 2  by  determine  the  major 
losses  and  minor  losses  through  the  control  volume  as  shown  in  Figure  2(a).  The  energy  equation  for  control 
volume  CV 2  is 


Pcb  -  Pace  =  {j^~KL2  ~  l)  , pVl  (7) 

where  Pcb  is  low  side  pressure  and  Pacc  is  accumulator  pressure.  Simplifying  Equation  (6),  (7)  and  substituting 
the  constitutive  equation  of  stack  yields 


P  = 


7  (  Pace  +  ~7 — 


X0  Ka 

CA, 


ch 


n  =  (y  (Kli  +  7 3KL2)  -  (7  -  1)  A2  -  1) 


(8) 


Equation  (8)  is  a  non-linear  first  order  differential  equation.  For  a  current  controlled  stack,  the  input  current 
i  of  the  stack  depends  on  the  characteristic  of  the  power  supply.  We  can  solve  the  differential  equation  of  the 
pressure  in  Equation  (8)  if  the  dynamics  of  the  input  current  are  known. 

Expansion  Stage 

Assuming  the  initial  pressure  is  P3  and  initial  displacement  £3,  the  equation  of  expansion  stage  becomes 

Vo  x 

P  =  P3~P  In^f - •  (9) 


Intake  Stage 

This  stage  is  similar  to  the  exhaustion  stage.  Figure  2(b)  shows  the  control  volume  CV3  in  dashed  lines. 
Using  the  same  procedure  as  in  the  exhaustion  stage,  the  flow  equations  in  the  intake  stage  are 


(10) 


where  0  =  L 3  +  l'j  and  Kl 3  is  the  loss  coefficient  in  the  control  volume. 

In  the  above  paragraph,  the  equations  of  flow  in  each  stage  are  derived  respectively.  Thus,  the  operations  of 
the  system  in  one  full  4-stage  cycle  can  be  modeled  by  combining  these  equations. 


2.1.1.  Summary  of  IVF  model 

In  the  IVF  model,  we  assume  a  viscous  incompressible  flow.  The  viscosity  is  taken  into  account  in  order  to 
predict  energy  loss  in  the  hydraulic  components.  Reynolds  number  is  computed  to  determine  the  flow  in  the 
system.  The  flow  equations  in  each  stage  is  derived  based  on  the  energy  equations.  By  assuming  the  fluid  is 
incompressible,  the  IVF  model  is  described  by  only  one  first  order  differential  equation  and  is  easy  to  solve  but 
it  might  overestimate  the  performance  by  ignoring  the  compressibility  of  the  fluid. 


2.2.  Compressible  viscous  flow  (CVF)  model 

In  this  model,  compressibility  of  the  fluid  will  be  modeled  in  the  full  cycle  of  the  four  stages.  Because  we  already 
assumed  compressibility  in  the  compression  stage  and  expansion  stage  in  the  development  of  the  IVF  model, 
equations  of  motion  of  IVF  model  in  these  two  stages  are  the  same  for  CVF  model.  Thus,  we  only  need  to  derive 
the  flow  equations  for  the  exhaustion  stage  and  intake  stage. 

The  system  is  considered  as  seven  components  of  fluid  as  shown  in  Figure  3.  To  maintain  the  simplicity  of  the 
model,  assumptions  and  approximations  are  made  based  on  observation  of  the  setup.  First  assumption  is  that 
we  assume  energy  losses,  including  major  losses  and  minor  losses  only  exist  in  the  pipes,  fittings  and  valves,  we 
assume  compressibility  of  fluid  only  occurs  in  pumping  chamber  and  high  side  of  the  cylinder,  which  corresponds 
to  part  1  and  3  respectively  in  Figure  3.  We  neglect  the  compressibility  in  component  4.  The  reason  is  that  it  is 
directly  connected  with  the  accumulator,  which  means  the  pressure  variation  is  small  and  thus  we  can  consider 
the  density  is  constant.  Because  the  cross  section  area  of  pipe  is  very  small  and  the  length  of  pipe  in  part  2  is 
made  short  as  possible  to  reduce  energy  losses,  we  can  assume  the  fluid  volume  in  the  pipe  is  relatively  small 


Pcb 


Figure  3:  Diagram  of  system  used  to  develop  CVF  model. 
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Table  2:  Key  equations  of  IVF  and  CVF  model  at  each  stage. 


compared  to  the  volume  in  the  pumping  chamber  and  hydraulic  cylinder.  Thus  we  neglect  the  mass  changes  of 
fluid  in  pipes  of  component  2. 

The  density  of  fluid  is  no  longer  a  constant  as  in  the  incompressible  model.  The  density-pressure  relationship 
of  fluid  can  be  derived  from  the  definition  of  bulk  modulus  in  Equation  (3)  as  following 


-  =  >• 

P  P 


(11) 


The  approach  of  this  modeling  is  to  incorporate  the  density  variation  into  the  incompressible  model.  The 
key  equations  of  the  hydraulic  components  in  each  stage  are  summarized  in  the  Table  2.  The  second  column  is 
the  component  corresponding  to  Figure  3.  The  energy  equations  of  CVF  model  is  the  same  as  IVF  model.  The 
only  difference  is  the  continuity  equation  in  the  pumping  chamber  and  hydraulic  cylinder,  where  density  change 
is  taken  into  account.  The  CVF  model  involves  coupled  differential  equations  and  can  be  numerically  solved. 


3.  EXPERIMENTAL  SETUP 

The  bench  top  experimental  setup  used  can  be  seen  in  Figure  1.  The  experimental  setup  incorporated  active 
valves  comprised  of  a  small  piezoceramic  stack  to  pump  fluid  and  move  a  hydraulic  cylinder.  Using  the  setup  as 
seen  in  Figure  1,  experiments  were  conducted  to  characterize  the  overall  actuator  performance.  The  piezoceramic 
stack  actuator  and  active  valves  were  all  driven  and  controlled  by  a  dSPACE  system  that  also  acquired  and 
recorded  the  data  from  the  sensors  used  in  the  experiments.  A  Sentec  pressure  sensor  is  used  to  measure 
the  pressure  inside  the  chamber.  A  Polytec  GmbH  laser  vibrometer  (model  OFV  303)  and  controller  (model 
OFV  3001)  were  used  to  measure  the  motion  of  a  target  connected  to  the  end  of  the  hydraulic  cylinder.  The 
laser  vibrometer  measured  the  cylinder  displacement  over  a  finite  amount  of  time  while  pulling  the  load. 


Simulation  of  velocity  of  load  at  50  Hz  frequency 


Simulation  of  output  power  at  at  50  Hz  frequency 
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(a)  Velocity  comparison  at  50  Hz 
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(b)  Power  comparison  at  50  Hz 
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Simulation  of  output  power  at  at  1 00  Hz  frequency 


(c)  Velocity  comparison  at  100  Hz 


(d)  Power  comparison  at  100  Hz 


Figure  4:  Comparison  of  simulated  and  experimental  velocity  and  power  for  a  single-ended  cylinder. 

4.  COMPARISON  OF  SIMULATIONS  AND  TEST  RESULTS 

4.1.  Simulation  parameters 

The  model  of  IVF  and  CVF  presented  derived  in  the  previous  section  will  be  used  to  predict  the  performance 
and  efficiency  calculations.  In  the  simulations  performed  in  this  paper  the  following  parameters  will  be  used: 

-  The  bulk  modulus  of  the  fluid  will  be  modeled  as  linear  function  of  pre-pressure  (10): 

0  =  48.8  P  +  15356,  /30  =  15356,  7  =  48.8  (12) 

This  linear  model  is  characterized  by  closed  chamber  tests  where  stiffness  of  the  fluid  is  measured. 

-  The  density  of  fluid  is  assumed  to  be  p  =  900 (kg/m3)  ,  which  is  a  typical  density  for  mineral  oil.  The 
viscosity  v  =  3  x  l(U5(m/s2). 

-  Two  setups  were  tested.  The  first  one  is  the  single-ended  cylinder  setup  powered  by  DSM  2  Level  PZT 
drive  whole  maximum  current  is  4.2  A.  The  other  is  the  double-ended  cylinder  powered  by  DSM  3-channel  power 
supply  whose  maximum  current  is  32.5  A. 

4.2.  Single-ended  cylinder  setup 

The  fluid  was  pre-pressured  at  230  Psi  for  the  first  set  of  tests.  Because  the  maximum  operation  frequency  of 
the  DSM  2  amplifier  is  100  Hz,  the  system  was  tested  at  frequencies  from  20  Hz  up  to  100  Hz.  Through  these 
tests,  a  timing  scheme  of  50%  stack  duty  cycle  and  40%  valve  duty  cycle  with  5%  outlet  valve  offset  and  55% 
inlet  valve  offset  is  adopted.  The  pressure  in  the  pumping  chamber,  the  displacement  of  load  and  displacement 
of  the  piston  are  measured. 

Figure  4  shows  the  comparison  of  simulation  results  and  experimental  data  for  50Hz  and  100  Hz  respectively. 
The  CVF  model  is  much  more  close  to  experimental  data  than  IVF  model  in  terms  of  both  velocity  and  power. 
For  loads  up  to  15  Kg,  the  simulation  results  of  CVF  model  consist  with  the  experimental  data  pretty  well.  It 
is  noted  the  simulation  velocity  of  both  model  are  very  like  a  straight  line.  But  the  slope  of  the  CVF  model  is 
steeper  than  the  slope  of  IVF.  Accordingly,  the  difference  of  the  two  models  becomes  more  obvious  when  the 
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Figure  5.  Comparison  of  simulated  and  experimental  output  power  and  velocity  under  130  V  stack  excitation  for  a 
double-ended  cylinder.  (a),(b),(c),(d):  Velocity  at  0,  6.84,  12.82,  17.45  Kg. 


load  is  larger  than  7  Kg.  The  reason  can  be  explained  by  the  compressibility  of  the  fluid.  When  a  larger  load  is 
attached  to  the  cylinder,  the  back  pressure  is  increased,  which  reduces  the  stroke  of  the  piezoelectric  actuator. 
Thus  the  displacement  of  the  load  also  decreases  accordingly.  Furthermore,  due  to  the  compressibility  of  the 
load,  additional  decrease  of  load  displacement  is  inevitable.  When  the  high  pressure  fluid  in  the  chamber  flows 
into  the  low  pressure  cylinder,  the  volume  will  increases  due  to  pressure  drop.  A  heavy  load  results  in  higher 
pressure  in  the  cylinder  and  thus  lower  pressure  drop,  the  effect  of  the  volume  increases  by  the  pressure  drop  is 
reduced.  This  additional  decrease  of  load  displacement  is  not  taken  into  account  for  the  IVF  models  because  they 
assume  an  incompressible  fluid.  Hence,  the  simulation  results  of  IVF  model  overestimates  system  performance 
under  large  loads.  It  is  also  noted  that  for  load  heavier  than  15Kg,  the  difference  between  the  CVF  model  and 
experimental  data  increases.  As  the  load  increases  over  15Kg,  the  measured  velocity  decreases  very  fast  for  both 
frequencies.  We  believe  this  fast  drop  of  performance  for  heavy  load  is  because  the  leakage  problem  of  the  valves. 

4.3.  Double-ended  cylinder  setup 

The  double-ended  cylinder  test  bench  was  built  in  order  to  achieve  better  performance  and  it  was  tested  using 
a  new  amplifier  that  enables  operation  in  the  kilo-hertz  range.  Different  loads  were  attached  with  the  output 
cylinder  when  operated  in  a  series  of  frequencies. 

The  simulation  of  velocity  versus  frequency  at  different  loads  are  plotted  in  the  following  figures.  Due  to 
the  hydraulic  amplification  mechanism  of  the  active  valves,  we  expect  some  delay  might  be  introduced.  In 
the  simulations,  a  0.4  ms  valve  transition  is  assumed.  Figure  5compares  respectively  the  velocities  and  power 
respectively  the  velocities  of  the  IVF  and  CVF  models  with  experimental  results  at  130  V  peak  voltages.  The 
CVF  model  consists  with  the  experimental  model  better  than  IVF  model.  The  performance  of  the  system  topped 
at  about  140  Hz.  Once  again,  we  notice  that  under  heavy  load,  the  large  difference  exists  between  the  model 
simulation  and  experimental  results.  We  also  attribute  this  discrepancy  to  the  leakage  problem  that  is  not 
accounted  in  the  model  simulations. 

We  observed  that  as  the  load  increases  the  IVF  model  overestimates  the  system  performance  as  it  does  in 
the  single-ended  cylinder  setup,  which  can  also  be  attributed  to  compressibility  of  fluid.  Another  observation  is 
that  the  difference  between  the  IVF  and  CVF  becomes  more  obvious  as  frequency  increases  in  each  plot.  The 


reason  of  this  discrepancy  is  that  CVF  model  assumes  a  linear  bulk  modulus  to  pressure  model  while  the  IVF 
model  assumes  an  infinite  modulus.  As  the  frequency  goes  higher,  because  the  fluid  can  not  flow  fast  enough  to 
fill  the  chamber  to  the  pre-pressure,  the  stiffness  of  the  fluid  decreased  due  to  reduce  of  bulk  modulus.  A  lower 
stiffness  of  fluid  results  in  lower  dynamic  pressure  and  smaller  stroke  of  the  load  in  one  cycle.  This  effect  ends 
up  to  a  big  difference  between  the  IVF  model  and  CVF  model  when  the  pump  is  operated  at  frequency  near  200 
Hz.  It  indicates  the  compressibility  of  fluid  is  important  to  the  performance  of  the  system. 

5.  SUMMARY  AND  CONCLUSIONS 

Two  steady  flow  models  are  presented  in  this  paper  based  on  different  approximation  of  the  fluid.  The  in¬ 
compressible  viscous  flow  (IVF)  model  is  developed  in  favor  of  minimizing  computation  effort  by  ignoring  the 
compressibility  of  the  hydraulic  fluid.  The  compressible  viscous  flow  (CVF)  model  is  then  presented  in  order  to 
achieve  more  accurate  results.  While  compressibility  of  fluid  is  taken  into  account,  assumptions  and  approxima¬ 
tions  are  carefully  made  to  maintain  the  simplicity  for  numerical  simulations. 

A  single-ended  cylinder  setup  and  a  double-ended  cylinder  setup  were  developed  and  tested.  The  test  results 
were  compared  with  theoretical  simulations.  The  comparison  revealed  the  IVF  model  deviated  from  the  CVF 
model  under  higher  loads,  which  indicates  the  compressibility  of  the  fluid  is  important  for  the  performance  of  the 
stack.  Increasing  pre-pressure  thus  increasing  bulk  modulus  of  fluid  can  potentially  improve  the  performance. 
Under  low  loads,  it  is  accurate  to  use  IVF  model  to  predict  performance.  The  CVF  model  provide  better  results 
under  higher  loads  and  the  CVF  model  consists  with  test  results  well  for  loads  up  to  20Kg.  But  when  the 
load  is  increased  above  that  limit,  the  performance  of  the  system  drops  more  sharply  than  the  model  predicts. 
Combined  with  the  reverse  movement  of  load  observed  under  large  loads  and  the  pressure  drop  in  closed  chamber 
figuration,  we  pretty  confidently  attribute  the  reason  to  the  leakage  problem  through  the  valves. 

The  simulation  results  and  experimental  tests  indicate  that  friction  in  the  hydraulic  circuits  is  a  major 
limitation  of  the  system  performance  when  going  to  higher  frequencies.  The  maximum  stroke  of  piston  decreases 
and  peak  pressure  drops,  resulting  in  a  reduction  of  performance.  The  system  performance  is  maximized  when 
operated  at  around  140  Hz.  Based  on  this  observation,  a  possible  effective  way  to  optimize  the  system  is  to 
reduce  friction  and  head  losses  in  the  hydraulic  circuits. 
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